Scalable, Cost-Effective, and Decentralized DNA Barcoding with Oxford Nanopore Sequencing.

DNA barcodes are useful in biodiversity research, but sequencing barcodes with dye termination methods ("Sanger sequencing") has been so time-consuming and expensive that DNA barcodes are not as widely used as they should be. Fortunately, MinION sequencers from Oxford Nanopore Technologies have recently emerged as a cost-effective and efficient alternative for barcoding. MinION barcodes are now suitable for large-scale species discovery and enable specimen identification when the target species are represented in barcode databases. With a MinION, it is possible to obtain 10,000 barcodes from a single flow cell at a cost of less than 0.10 USD per specimen. Additionally, a Flongle flow cell can be used for small projects requiring up to 300 barcodes (0.50 USD per specimen). We here describe a cost-effective laboratory workflow for obtaining tagged amplicons, preparing ONT libraries, sequencing amplicon pools, and analyzing the MinION reads with the software ONTbarcoder. This workflow has been shown to yield highly accurate barcodes that are 99.99% identical to Sanger barcodes. Overall, we propose that the use of MinION for DNA barcoding is an attractive option for all researchers in need of a cost-effective and efficient solution for large-scale species discovery and specimen identification.

De novo sequencing, diploid assembly, and annotation of the black carpenter ant, Camponotus pennsylvanicus, and its symbionts by one person for $1000, using nanopore sequencing.

The black carpenter ant (Camponotus pennsylvanicus) is a pest species found widely throughout North America. From a single individual I used long-read nanopore sequencing to assemble a phased diploid genome of 306 Mb and 60X coverage, with quality assessed by a 97.0% BUSCO score, improving upon other ant assemblies. The mitochondrial genome reveals minor rearrangements from other ants. The reads also allowed assembly of parasitic and symbiont genomes. I include a complete Wolbachia bacterial assembly with a size of 1.2 Mb, as well as a commensal symbiont Blochmannia pennsylvanicus, at 791 kb. DNA methylation and hydroxymethylation were measured at base-pair resolution level from the same reads and confirmed extremely low levels seen in the Formicidae family. There was moderate heterozygosity, with 0.16% of bases being biallelic from the parental haplotypes. Protein prediction yielded 14 415 amino acid sequences with 95.8% BUSCO score and 86% matching to previously known proteins. All assemblies were derived from a single MinION flow cell generating 20 Gb of sequence for a cost of $1047 including consumable reagents. Adding fixed costs for equipment brings the total for an ant-sized genome to less than $5000. All analyses were performed in 1 week on a single desktop computer.

Creating reference animal genomes is typically a large, expensive process. Here I sequenced the genome of the black carpenter ant for only $1000 as a sole researcher in just one week. Along with the nuclear genome, I assembled the mitochondrial genome and two commensal bacteria species living within the ant. Nanopore technology also enabled epigenetic measurements from the same ant and replicated other studies showing very low DNA methylation. The reference genome compared favorably to other ant species in continuity and protein prediction accuracy. This method will allow other low-resource labs to create high quality genome assemblies with a low cost.

Hackflex: low-cost, high-throughput, Illumina Nextera Flex library construction.

We developed a low-cost method for the production of Illumina-compatible sequencing libraries that allows up to 14 times more libraries for high-throughput Illumina sequencing to be generated for the same cost. We call this new method Hackflex. The quality of library preparation was tested by constructing libraries from Escherichia coli MG1655 genomic DNA using either Hackflex, standard Nextera Flex (recently renamed as Illumina DNA Prep) or a variation of standard Nextera Flex in which the bead-linked transposase is diluted prior to use. In order to test the library quality for genomes with a higher and a lower G+C content, library construction methods were also tested on Pseudomonas aeruginosa PAO1 and Staphylococcus aureus ATCC 25923, respectively. We demonstrated that Hackflex can produce high-quality libraries and yields a highly uniform coverage, equivalent to the standard Nextera Flex kit. We show that strongly size-selected libraries produce sufficient yield and complexity to support de novo microbial genome assembly, and that assemblies of the large-insert libraries can be much more contiguous than standard libraries without strong size selection. We introduce a new set of sample barcodes that are distinct from standard Illumina barcodes, enabling Hackflex samples to be multiplexed with samples barcoded using standard Illumina kits. Using Hackflex, we were able to achieve a per-sample reagent cost for library prep of A$7.22 (Australian dollars) (US $5.60; UK £3.87, £1=A$1.87), which is 9.87 times lower than the standard Nextera Flex protocol at advertised retail price. An additional simple modification and further simplification of the protocol by omitting the wash step enables a further price reduction to reach an overall 14-fold cost saving. This method will allow researchers to construct more libraries within a given budget, thereby yielding more data and facilitating research programmes where sequencing large numbers of libraries is beneficial.

Screening of CNVs using NGS data improves mutation detection yield and decreases costs in genetic testing for hereditary cancer.

INTRODUCTION: Germline CNVs are important contributors to hereditary cancer. In genetic diagnostics, multiplex ligation-dependent probe amplification (MLPA) is commonly used to identify them. However, MLPA is time-consuming and expensive if applied to many genes, hence many routine laboratories test only a subset of genes of interest. METHODS AND RESULTS: We evaluated a next-generation sequencing (NGS)-based CNV detection tool (DECoN) as first-tier screening to decrease costs and turnaround time and expand CNV analysis to all genes of clinical interest in our diagnostics routine. We used DECoN in a retrospective cohort of 1860 patients where a limited number of genes were previously analysed by MLPA, and in a prospective cohort of 2041 patients, without MLPA analysis. In the retrospective cohort, 6 new CNVs were identified and confirmed by MLPA. In the prospective cohort, 19 CNVs were identified and confirmed by MLPA, 8 of these would have been lost in our previous MLPA-restricted detection strategy. Also, the number of genes tested by MLPA across all samples decreased by 93.0% in the prospective cohort. CONCLUSION: Including an in silico germline NGS CNV detection tool improved our genetic diagnostics strategy in hereditary cancer, both increasing the number of CNVs detected and reducing turnaround time and costs.

Evaluation of a high-throughput, cost-effective Illumina library preparation kit.

Library preparation for high-throughput sequencing applications is a critical step in producing representative, unbiased sequencing data. The iGenomX Riptide High Throughput Rapid Library Prep Kit purports to provide high-quality sequencing data with lower costs compared to other Illumina library kits. To test these claims, we compared sequence data quality of Riptide libraries to libraries constructed with KAPA Hyper and NEBNext Ultra. Across several single-source genome samples, mapping performance and de novo assembly of Riptide libraries were similar to conventional libraries prepared with the same DNA. Poor performance of some libraries resulted in low sequencing depth. In particular, degraded DNA samples may be challenging to sequence with Riptide. There was little cross-well plate contamination with the overwhelming majority of reads belong to the proper source genomes. The sequencing of metagenome samples using different Riptide primer sets resulted in variable taxonomic assignment of reads. Increased adoption of the Riptide kit will decrease library preparation costs. However, this method might not be suitable for degraded DNA.

Endonuclease enrichment TAPS for cost-effective genome-wide base-resolution DNA methylation detection.

Whole genome base-resolution methylome sequencing allows for the most comprehensive analysis of DNA methylation, however, the considerable sequencing cost often limits its applications. While reduced representation sequencing can be an affordable alternative, over 80% of CpGs in the genome are not covered. Building on our recently developed TET-assisted pyridine borane sequencing (TAPS) method, we here described endonuclease enrichment TAPS (eeTAPS), which utilizes dihydrouracil (DHU)-cleaving endonuclease digestion of TAPS-converted DNA to enrich methylated CpG sites (mCpGs). eeTAPS can accurately detect 87% of mCpGs in the mouse genome with a sequencing depth equivalent to 4× whole genome sequencing. In comparison, reduced representation TAPS (rrTAPS) detected less than 4% of mCpGs with 2.5× sequencing depth. Our results demonstrate eeTAPS to be a new strategy for cost-effective genome-wide methylation analysis at single-CpG resolution that can fill the gap between whole-genome and reduced representation sequencing.

The profiling of microbiota in vaginal swab samples using 16S rRNA gene sequencing and IS-pro analysis.

BACKGROUND: 16S rRNA gene sequencing is currently the most common way of determining the composition of microbiota. This technique has enabled many new discoveries to be made regarding the relevance of microbiota to the health and disease of the host. However, compared to other diagnostic techniques, 16S rRNA gene sequencing is fairly costly and labor intensive, leaving room for other techniques to improve on these aspects. RESULTS: The current study aimed to compare the output of 16S rRNA gene sequencing to the output of the quick IS-pro analysis, using vaginal swab samples from 297 women of reproductive age. 16S rRNA gene sequencing and IS-pro analyses yielded very similar vaginal microbiome profiles, with a median Pearson's R2 of 0.97, indicating a high level of similarity between both techniques. CONCLUSIONS: We conclude that the results of 16S rRNA gene sequencing and IS-pro are highly comparable and that both can be used to accurately determine the vaginal microbiota composition, with the IS-pro analysis having the benefit of rapidity.

Resizing reaction volumes for the ForenSeq™ DNA Signature Prep kit library preparation.

The introduction of next generation sequencing (NGS; also known as massively parallel sequencing) technology in the field of forensic genetics has been welcomed by the scientific community, above all because it complements the weaknesses of capillary electrophoresis (CE) in the analysis of genetic markers, such as single nucleotide polymorphism (SNP) typing. However, one of the main obstacles to its adoption does not seem to be the cost of the instrumentation, but rather the cost of the NGS library preparation kits. With the aim of reducing the cost of library preparation without compromising the quality of the results, we tried to scale down reaction volumes for the first two polymerase chain reactions in the amplification and enrichment phases of the targeted loci of library preparation using the ForenSeq™ DNA Signature Prep kit. We used 1 µL templated DNA input to a concentration of 1 ng/µL, instead of the 5 µL at 0.2 ng/µL recommended by the manufacturer. Our findings indicate that reduction of the library preparation volume using the ForenSeq™ DNA Signature Prep kit did not interfere with the quality and reproducibility of the DNA profiles obtained and can help lower the overall cost of NGS.

A Novel Next-Generation Sequencing-Based Approach for Concurrent Detection of Mitochondrial DNA Copy Number and Mutation.

Numerous studies have identified essential contributions of altered mitochondrial DNA (mtDNA) copy number and mutations in many common disorders, including cancer. To date, capture-based next-generation sequencing (NGS) has been widely applied to detect mtDNA mutations, although it lacks the ability to assess mtDNA copy number. The current strategy for quantifying mtDNA copy number relies mainly on real-time quantitative PCR, which is limited in degraded samples. A novel capture-based NGS approach was developed using both mtDNA and nuclear DNA probes to capture target fragments, enabling simultaneous detection of mtDNA mutations and copy number in different sample types. First, the impact of selecting reference genes on mtDNA copy number calculation was evaluated, and finally, 3 nuclear DNA fragments of 4000 bp were selected as an internal reference for detection. Then, the effective application of this approach was verified in DNA samples of formalin-fixed, paraffin-embedded specimens and body fluids, indicating the widespread applicability. This approach showed more accurate and stable results in detecting mtDNA copy number compared with real-time quantitative PCR in degraded DNA samples. Moreover, data indicated this approach had good reproducibility in detecting both mtDNA copy number and mutations among three sample types. Altogether, a versatile and cost-effective capture-based NGS approach has been developed for concurrent detection of mtDNA copy number and mutations, which has numerous applications in research and diagnosis.

Cost-effectiveness of the CFTR gene-sequencing test for asymptomatic carriers in the Colombian population / [Costo-efectividad de la prueba de secuenciación del gen CFTR para portadores asintomáticos en la población colombiana].

Introduction: Cystic fibrosis is an autosomal recessive genetic disease classified as a highcost orphan disease. Objective: To determine the cost-effectiveness ratio of the diagnostic test for the CFTR gene-sequencing in asymptomatic family carriers in the first, second, and third degree of consanguinity. Materials and methods: We conducted a systematic search evaluating operative characteristics of the diagnostic test and decision-tree models in cost-effectiveness studies. A decision-tree model was elaborated taking prevention of future conceptions as a unit of analysis. We obtained the costs of the disease from the high-cost report of the Ministerio de Salud y Protección Social. The costs of the test were referenced by national laboratories. We carried out a deterministic and probabilistic sensitivity analysis with a third-payer perspective and a one-year horizon. Results: An ICER of USD$ 5051.10 was obtained as the incremental cost for obtaining 10.89% more probability of avoiding the birth of a child with cystic fibrosis per screened couple. For family members in second and third degrees, the ICER was USD$ 19,380.94 and USD$ 55,913.53, respectively, evidenced when applying the GDP per capita. This technology was cost-effective in 39%, 61.18%, and 74.36% for 1, 2, and 3 GDP per capita in first degree of consanguinity relatives. Conclusions: The genetic test for the detection of CFTR gene carriers was cost-effective depending on the threshold of availability to pay and the assumptions and limitations established in the model.

Introducción. La fibrosis quística es una enfermedad genética de carácter autosómico recesivo clasificada como enfermedad huérfana de alto costo. Objetivo. Determinar la razón de costo-efectividad de la prueba diagnóstica de secuenciación del gen CFTR para los portadores asintomáticos familiares en primer, segundo y tercer grados de consanguinidad. Materiales y métodos. Se hizo una búsqueda sistemática sobre la evaluación de las características operativas de la prueba diagnóstica y los modelos de árbol de decisiones en estudios de costo-efectividad. Se elaboró un modelo de árbol de decisiones tomando como unidad de análisis la prevención de futuras concepciones. Los costos de la enfermedad se obtuvieron del reporte de alto costo del Ministerio de Salud de Colombia. Los costos de la prueba se obtuvieron de laboratorios nacionales. Se hizo un análisis de sensibilidad, determinístico y probabilístico, con la perspectiva del tercer pagador y horizonte a un año. Resultados. Se obtuvo una razón incremental de costo-efectividad (RICE) de USD$5.051,10 por obtener 10,89 % más de probabilidades de evitar el nacimiento de un niño enfermo con fibrosis quística por pareja. Para los familiares de segundo y tercer grados, se encontró una RICE de USD$ 19.380,94 y USD$ 55.913,53, respectivamente, al aplicar el PIB per cápita. Esta tecnología fue costo-efectiva en 39 %, 61,18 % y 74,36 % para 1, 2 y 3 PIB per cápita en familiares de primer grado de consanguinidad. Conclusiones. La prueba genética de detección de portadores del gen CFTR resultó costo-efectiva dependiendo del umbral de la disponibilidad de pagar, y de los supuestos y limitaciones establecidas en el modelo.

A cost-effective approach to DNA methylation detection by Methyl Sensitive DArT sequencing.

Several studies suggest the relation of DNA methylation to diseases in humans and important phenotypes in plants drawing attention to this epigenetic mark as an important source of variability. In the last decades, several methodologies were developed to assess the methylation state of a genome. However, there is still a lack of affordable and precise methods for genome wide analysis in large sample size studies. Methyl sensitive double digestion MS-DArT sequencing method emerges as a promising alternative for methylation profiling. We developed a computational pipeline for the identification of DNA methylation using MS-DArT-seq data and carried out a pilot study using the Eucalyptus grandis tree sequenced for the species reference genome. Using a statistic framework as in differential expression analysis, 72,515 genomic sites were investigated and 5,846 methylated sites identified, several tissue specific, distributed along the species 11 chromosomes. We highlight a bias towards identification of DNA methylation in genic regions and the identification of 2,783 genes and 842 transposons containing methylated sites. Comparison with WGBS, DNA sequencing after treatment with bisulfite, data demonstrated a precision rate higher than 95% for our approach. The availability of a reference genome is useful for determining the genomic context of methylated sites but not imperative, making this approach suitable for any species. Our approach provides a cost effective, broad and reliable examination of DNA methylation profile on MspI/HpaII restriction sites, is fully reproducible and the source code is available on GitHub (https://github.com/wendelljpereira/ms-dart-seq).

Targeted-Sequencing Workflows for Comprehensive Drug Resistance Profiling of Mycobacterium tuberculosis Cultures Using Two Commercial Sequencing Platforms: Comparison of Analytical and Diagnostic Performance, Turnaround Time, and Cost.

BACKGROUND: The emergence of Mycobacterium tuberculosis with complex drug resistance profiles necessitates a rapid and comprehensive drug susceptibility test for guidance of patient treatment. We developed two targeted-sequencing workflows based on Illumina MiSeq and Nanopore MinION for the prediction of drug resistance in M. tuberculosis toward 12 antibiotics. METHODS: A total of 163 M. tuberculosis isolates collected from Hong Kong and Ethiopia were subjected to a multiplex PCR for simultaneous amplification of 19 drug resistance-associated genetic regions. The amplicons were then barcoded and sequenced in parallel on MiSeq and MinION in respective batch sizes of 24 and 12 samples. A web-based bioinformatics pipeline, BacterioChek-TB, was developed to translate the raw datasets into clinician-friendly reports. RESULTS: Both platforms successfully sequenced all samples with mean read depths of 1,127× and 1,649×, respectively. The variant calling by MiSeq and MinION could achieve 100% agreement if variants with an allele frequency of <40% reported by MinION were excluded. Both workflows achieved a mean clinical sensitivity of 94.8% and clinical specificity of 98.0% when compared with phenotypic drug susceptibility test (pDST). Turnaround times for the MiSeq and MinION workflows were 38 and 15 h, facilitating the delivery of treatment guidance at least 17-18 days earlier than pDST, respectively. The higher cost per sample on the MinION platform ($71.56) versus the MiSeq platform ($67.83) was attributed to differences in batching capabilities. CONCLUSION: Our study demonstrates the interchangeability of MiSeq and MinION platforms for generation of accurate and actionable results for the treatment of tuberculosis.

Cost-effectiveness analysis of pretreatment screening for NUDT15 defective alleles.

BACKGROUND: Nucleotide triphosphate diphosphatase (NUDT15) genetic testing in addition to thiopurine methyl transferase (TPMT) is recommended to reduce the incidence of adverse severe myelotoxicity episodes induced by thiopurines. OBJECTIVE: We assessed the cost-effectiveness ratio of combined screening for TMPT and NUDT15 defective alleles by genotyping or next-generation sequencing (NGS) using TPMT genotyping as the reference. Because of the genetic differences in thiopurine toxicity, we tested the screening strategies on individuals of Caucasian and Asian descent. METHODS: A decision tree compared conventional TPMT genotyping with combined TPMT/NUDT15 genotyping or NGS using a Monte-Carlo microsimulation model of patients with inflammatory bowel disease. The main outcome was the incremental cost-effectiveness ratios (ICER) with effectiveness being one averted severe myelotoxicity requiring hospitalization. RESULTS: The mean estimated cost of the TPMT genotyping for one year is twice in Asian compared with Caucasian patients (980 euro/patient versus 488 euro/patient), and the effectiveness of TPMT genotyping in Caucasian avoided 43 severe myelosuppressions per 10 000 patients over a year compared with 3.6 per 10 000 patients in Asian. Combined TPMT/NUDT15 genotyping compared with TPMT genotyping had an ICER of 7 491 281 euro per severe myelotoxicity averted in Caucasian, compared to 619 euro in Asian. The ICER of the NGS-based screening strategy is disproportionally high compared with genotyping, irrespective of ethnic descent. CONCLUSION: With a low cost-effectiveness threshold, combined screening for NUDT15 and TPMT defective alleles is cost-effective compared to TMPT screening alone in patients of Asian descent, but is unrealistic from a cost-effectiveness point of view in Caucasians.

The cost trajectory of the diagnostic care pathway for children with suspected genetic disorders.

PURPOSE: This study describes the cost trajectory of the standard diagnostic care pathway for children with suspected genetic disorders in British Columbia, Canada. METHODS: Average annual per-patient costs were estimated using medical records review and a caregiver survey for a cohort of 498 children referred to BC Children's and Women's Hospitals (C&W) with unexplained intellectual disability (the TIDE-BC study) and families enrolled in the CAUSES study, which offered diagnostic genome-wide sequencing (GWS; exome and genome sequencing) to 500 families of children with suspected genetic disorders. RESULTS: Direct costs peaked in the first year of patients' diagnostic odyssey, with an average of C$2257 per patient (95% confidence interval [CI] C$2074, C$2441) for diagnostic testing and C$631 (95% CI C$543, C$727) for specialist consultations at C&W. In subsequent years, direct costs accrued at a constant rate, with an estimated annual per-patient cost of C$511 (95% CI C$473, C$551) for diagnostic testing and C$334 (95% CI C$295, C$369) for consultations at C&W. Travel costs and caregiver productivity loss associated with attending diagnosis-related physician appointments averaged C$1907/family/year. CONCLUSIONS: The continuing long-term accrual of costs by undiagnosed patients suggests that economic evaluations of diagnostic GWS services should use longer time horizons than have typically been used.

DNA enrichment and tagmentation method for species-level identification and strain-level differentiation using ON-rep-seq.

Despite the massive developments within culture-independent methods for detection of microorganisms during the last decade, culture-based methods remain a cornerstone in microbiology. Yet, the problem of rapid, accurate and inexpensive identification of bacterial isolates down to species/strain level remains unresolved. We have developed a new method for bacterial DNA enrichment and tagmentation allowing fast (<24 h) and cost-effective species level identification and strain level differentiation using the MinION portable sequencing platform (ON-rep-seq). DNA library preparation for 96 isolates takes less than 5 h and ensures highly reproducible distribution of reads that can be used to generate strain level specific read length counts profiles (LCp). We have developed a pipeline that by correcting reads error within peaks of LCp generates a set of high quality (>99%) consensus reads. Whereas, the information from high quality reads is used to retrieve species level taxonomy, comparison of LCp allows for strain level differentiation.

Comparison of long-read sequencing technologies in the hybrid assembly of complex bacterial genomes.

Illumina sequencing allows rapid, cheap and accurate whole genome bacterial analyses, but short reads (<300 bp) do not usually enable complete genome assembly. Long-read sequencing greatly assists with resolving complex bacterial genomes, particularly when combined with short-read Illumina data (hybrid assembly). However, it is not clear how different long-read sequencing methods affect hybrid assembly accuracy. Relative automation of the assembly process is also crucial to facilitating high-throughput complete bacterial genome reconstruction, avoiding multiple bespoke filtering and data manipulation steps. In this study, we compared hybrid assemblies for 20 bacterial isolates, including two reference strains, using Illumina sequencing and long reads from either Oxford Nanopore Technologies (ONT) or SMRT Pacific Biosciences (PacBio) sequencing platforms. We chose isolates from the family Enterobacteriaceae, as these frequently have highly plastic, repetitive genetic structures, and complete genome reconstruction for these species is relevant for a precise understanding of the epidemiology of antimicrobial resistance. We de novo assembled genomes using the hybrid assembler Unicycler and compared different read processing strategies, as well as comparing to long-read-only assembly with Flye followed by short-read polishing with Pilon. Hybrid assembly with either PacBio or ONT reads facilitated high-quality genome reconstruction, and was superior to the long-read assembly and polishing approach evaluated with respect to accuracy and completeness. Combining ONT and Illumina reads fully resolved most genomes without additional manual steps, and at a lower consumables cost per isolate in our setting. Automated hybrid assembly is a powerful tool for complete and accurate bacterial genome assembly.

DiscoverY: a classifier for identifying Y chromosome sequences in male assemblies.

BACKGROUND: Although the Y chromosome plays an important role in male sex determination and fertility, it is currently understudied due to its haploid and repetitive nature. Methods to isolate Y-specific contigs from a whole-genome assembly broadly fall into two categories. The first involves retrieving Y-contigs using proportion sharing with a female, but such a strategy is prone to false positives in the absence of a high-quality, complete female reference. A second strategy uses the ratio of depth of coverage from male and female reads to select Y-contigs, but such a method requires high-depth sequencing of a female and cannot utilize existing female references. RESULTS: We develop a k-mer based method called DiscoverY, which combines proportion sharing with female with depth of coverage from male reads to classify contigs as Y-chromosomal. We evaluate the performance of DiscoverY on human and gorilla genomes, across different sequencing platforms including Illumina, 10X, and PacBio. In the cases where the male and female data are of high quality, DiscoverY has a high precision and recall and outperforms existing methods. For cases when a high quality female reference is not available, we quantify the effect of using draft reference or even just raw sequencing reads from a female. CONCLUSION: DiscoverY is an effective method to isolate Y-specific contigs from a whole-genome assembly. However, regions homologous to the X chromosome remain difficult to detect.